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Abstract. Many 2D lattice models of physical phenomena are conjectured to have con- 
formally invariant scaling limits: percolation, Ising model, self-avoiding polymers, . . . This 
has led to numerous exact (but non-rigorous) predictions of their scaling exponents and 
dimensions. We will discuss how to prove the conformal invariance conjectures, especially 
in relation to Schramm-Loewner Evolution. 
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1. Introduction 

For several 2D lattice models physicists were able to make a number of spectacu- 
lar predictions (non-rigorous, but very convincing) about exact values of various 
scaling exponents and dimensions. Many methods were employed (Coulomb Gas, 
Conformal Field Theory, Quantum Gravity) with one underlying idea: that the 
model at criticality has a continuum scaling limit (as mesh of the lattice goes to 
zero) and the latter is conformally invariant. Moreover, it is expected that there 
is only a one-parameter family of possible conformally invariant scaling limits, so 
universality follows: if the same model on different lattices (and sometimes at dif- 
ferent temperatures) has a conformally invariant scaling limit, it is necessarily the 
same. Indeed, the two limits belong to the same one-parameter family, and usually 
it directly follows that the corresponding parameter values coincide. 

Recently mathematicians were able to offer different, perhaps better, and cer- 
tainly more rigorous understanding of those predictions, in many cases providing 
proofs. The point which is perhaps still less understood both from mathematics 
and physics points of view is why there exists a universal conformally invariant 
scaling limit. However such behavior is supposed to be typical in 2D models at 
criticality: Ising, percolation, self-avoiding polymers; with universal conformally 
invariant curves arising as scaling limits of the interfaces. 

Until recently this was established only for the scaling limit of the 2D random 
walk, the 2D Brownian motion. This case is easier and somehow exceptional be- 
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cause of the Markov property. Indeed, Brownian motion was originally constructed 
by Wiener [46j , and its conformal invariance (which holds in dimension 2 only) was 
shown by Paul Levy [21] without appealing to random walk. Note also that un- 
like interfaces (which are often simple, or at most "touch" themselves), Brownian 
trajectory has many "transversal" self-intersections. 

For other lattice models even a rigorous formulation of conformal invariance 
conjecture seemed elusive. Considering percolation (a model where vertices of a 
graph are declared open independently with equal probability p - see the discus- 
sion below) at criticality as an example, Robert Langlands, Philippe Pouliot and 
Yvan Saint- Aubin in |20j studied numerically crossing probabilities (of events that 
there is an open crossing of a given rectangular shape). Based on experiments 
they concluded that crossing probabilities should have a universal (independent of 
lattice) scaling limit, which is conformally invariant (a conjecture they attributed 
to Michael Aizenman). Thus the limit of crossing probability for a rectangular do- 
main should depend on its conformal modulus only. Moreover an exact formula ([5]) 
using hypergeometric function was proposed by John Cardy in [9] based on Confor- 
mal Field Theory arguments. Later Lennart Carleson found that the formula has 
a particularly nice form for equilateral triangles, see [38j . These developments got 
many researchers interested in the subject and stimulated much of the subsequent 
progress. 

Rick Kenyon [16l [17] established conformal invariance of many observables 
related to dimer models (domino tilings), in particular to uniform spanning tree 
and loop erased random walk, but stopped short of constructing the limiting curves. 

In [32] , Oded Schramm suggested to study the scaling limit of a single interface 
and classified all possible curves which can occur as conformally invariant scal- 
ing limits. Those turned out to be a universal one-parameter family of SLE (k) 
curves, which are now called Schramm- Loewner Evolutions. The word "evolution" 
is used since the curves are constructed dynamically, by running classical Loewner 
Evolution with Brownian motion as a driving term. We will discuss one possible 
setup, chordal SLE (k) with parameter k € [0, oo), which provides for each simply- 
connected domain and boundary points a, b a measure n on curves from a to b 
inside tt. The measures a, b) are conformally invariant, in particular they are 
all images of one measure on a reference domain, say a half-plane C+. An exact 
definition appears below. 

In [3Tl EE] the conformal invariance was established for critical percolation on 
triangular lattice. Conformally invariant limit of the interface was identified with 
SLE (6), though its construction does not use SLE machinery. See also Federico 
Camia and Charles Newman's paper (8] for the details on subsequent construction 
of the full scaling limit. 

In [23] Greg Lawler, Oded Schramm and Wendelin Werner have shown that a 
perimeter curve of the uniform spanning tree converges to SLE (8) (and the related 
loop erased random walk - to SLE (2)) on a general class of lattices. Unlike the 
proof for percolation, theirs utilizes SLE in a substantial way. In [35] Oded 
Schramm and Scott Sheffield introduced a new model, Harmonic Explorer, where 
properties needed for convergence to SLE (4) are built in. 
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Despite the results for percolation and uniform spanning tree, the problem 
remained open for all other classical (spin and random cluster) 2D models, includ- 
ing percolation on other lattices. This was surprising given the abundance of the 
physics literature on conformal invariance. Perhaps most surprising was that the 
problem of a conformally invariant scaling limit remained open for the Ising model, 
since for the latter there are many exact and often rigorous results - see the books 

USE]. 

Recently we were able to work out the Ising case [39, 4"0[ I4T1 142] : 

Theorem 1. As lattice step goes to zero, interfaces in Ising and Ising random 
cluster models on the square lattice at critical temperature converge to SLE (3) and 
SLE(16/3) correspondingly. 

Computer simulations of these interfaces (Figures [2 [4]) as well as the defini- 
tion of the Ising models can be found below. Similarly to mentioned experiments 
for percolation, Robert Langlands, Marc Andre Lewis and Yvan Saint- Aubin con- 
ducted in [21] numerical studies of crossing probabilities for the Ising model at 
critical temperature. A modification of the theorem above relating interfaces to 
SLE's (with drifts) in domains with five marked boundary points allows a rigorous 
setup for establishing their conjectures. 

The proof is based on showing that a certain Fermionic lattice observable (or 
rather two similar ones for spin and random cluster models) is discrete analytic and 
solves a particular covariant Riemann Boundary Value Problem. Hence its limit is 
conformally covariant and can be calculated exactly. The statement is interesting 
in its own right, and can be used to study spin correlations. The observable studied 
has more manifest physics meaning than one in our percolation paper [38j . 

The methods lead to some progress in fairly general families of random cluster 
and 0(n) models, and not just on square lattices. In particular, besides Ising cases, 
they seem to suggest new proofs for all other known cases (i.e. site percolation on 
triangular lattice and uniform spanning tree). 

In this note we will discuss this proof and general approach to scaling limits 
and conformal invariance of interfaces in the SLE context. We will also state some 
of the open questions and speculate on how one should approach other models. 

We omit many aspects of this rich subject. We do not discuss the general 
mathematical theory of SLE curves or their connections to physics, for which 
interested reader can consult expository works [10l |T5l [44] and a book [22] . We 
do not mention the question of how to deduce the values of scaling exponents 
for lattice models with SLE help once convergence is known. It was explored in 
some detail only for percolation [43j , where convergence is known and the required 
(difficult) estimates were already in place thanks to Harry Kesten [19"] . We also 
restrict ourselves to one interface, whereas one can study the collection of all loops 
(cf. exposition [45]), and many of our considerations transfer to the loop soup 
observables. Finally, there are many other open questions related to conformal 
invariance, some of which are discussed in Oded Schramm's paper [34] in these 
proceedings. 
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2. Lattice models 

We focus on two families of lattice models which have nice "loop representations" . 
Those families include or are closely related to most of the "important" models, 
including percolation, Ising, Potts, spherical (or O(n)), Fortuin-Kasteleyn (or ran- 
dom cluster), self-avoiding random walk, and uniform spanning tree models. For 
their interrelations and for the discussion of many other relevant models one can 
consult the books [51 H21 HS1 [23 • We also omit many references which can be found 
there. 

There are various ways to understand the existence of the scaling limit and its 
conformal invariance. One can ask for the full picture, which can be represented 
as a loop collection (representing all cluster interfaces), random height function 
(changing by ±1 whenever we cross a loop), or some other object. It however 
seems desirable to start with a simpler problem. 

One can start with observables (like correlation functions, crossing probabili- 
ties), for which it is easier to make sense of the limit: there should exist a limit of a 
number sequence which is a conformal invariant. Though a priori it might seem to 
be a weaker goal than constructing a full scaling limit, there are indications that 
to obtain the full result it might be sufficient to analyze just one observable. 

We will discuss an intermediate goal to analyze the law of just one interface, 
explain why working out just one observable would be sufficient, and give details 
on how to find an observable with a conformally invariant limit. To single out 
one interface, we consider a model on a simply connected domain with Dobrushin 
boundary conditions (which besides many loop interfaces enforce existence of an 
interface joining two boundary points a and b). We omit the discussion of the full 
scaling limit, as well as models on Riemann surfaces and with different boundary 
conditions. 

2.1. Percolation. Perhaps the simplest model (to state) is Bernoulli perco- 
lation on the triangular lattice. Vertices are declared open or closed (grey or white 
in Figure [1]) independently with probabilities p and (1 — p) correspondingly. The 
critical value is p — p c — 1/2 - see [18l E] , in which case all colorings are equally 
probable. 

Then each configuration can be represented by a collection of interfaces - loops 
which go along the edges of the dual hexagonal lattice and separate open and 
closed vertices. 

We want to distinguish one particular interface, and to this effect we introduce 
Dobrushin boundary conditions: we take two boundary points a and b in a simply 
connected fi (or rather its lattice approximation), asking the counterclockwise arc 
ab to be grey and the counterclockwise arc ba to be white. This enforces existence 
of a single non-loop interface which runs from a to b. The "loop gas" formulation 
of our model is that we consider all collections of disjoint loops plus a curve from 
a to b on hexagonal lattice with equal probability. 

For each value of the lattice step e > we approximate a given domain O by 
a lattice domain, which leads to a random interface, that is a probability measure 
/i e on curves (broken lines) running from a to b. The question is whether there is a 
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Figure 1. Critical site percolation on triangular lattice superimposed over a rectangle. 
Every site is grey or white independently with equal probability 1/2. Dobrushin boundary 
conditions (grey on lower and left sides, white on upper and right sides) produce an 
interface from the upper left corner a to the lower right corner b. The law of the interface 
converges to SLE (6) when lattice step goes to zero, while the rectangle is fixed. 

limit measure /i = fJ,(Q, a, b) on curves and whether it is conformally invariant. To 
make sense of the limit we consider the curves with uniform topology generated by 
parameterizations (with distance between 71 and 72 being inf ||/i — /2II00 where the 
infimum is taken over all parameterizations /1, of 71, 72), and ask for weak-* 
convergence of the measures /j. e . 

2.2. O(n) and loop models. Percolation turns out to be a particular case 
of the loop gas model which is closely related (via high-temperature expansion) to 
0(n) (spherical) model. We consider configurations of non- intersecting simple 
loops and a curve running from a to 6 on hexagonal lattice inside domain as 
for percolation in Figure [T] But instead of asking all configurations to be equally 
likely, we introduce two parameters: loop-weight n > and edge-weight x > 0, 
and ask that probability of a configuration is proportional to 

^# loops length of loops 

The vertices not visited by loops are called monomers. Instead of weighting edges 
by x one can equivalently weight monomers by l/x. 

We are interested in the range n £ [0, 2] (after certain modifications n £ [—2, 2] 
would work), where conformal invariance is expected (other values of n have differ- 
ent behavior). It turns out that there is a critical value x c (n), such that the model 
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exhibits one critical behavior at x c (n) and another on the interval (x c (n), +00), 
corresponding to "dilute" and "dense" phases (when in the limit the loops are 
simple and non-simple correspondingly). 

Bernard Nienhuis [2H1 US] proposed the following conjecture, supported by 
physics arguments: 

Conjecture 2. The critical value is given by 

1 



x c {n) 



Note that though for all x E (x c (n),oo) the critical behavior (and the scaling limit) 



are conjecturally the same, the related value x c (n) = l/y 2 — \/2 — n turns out to 
be distinguished in some ways. 

The criticality was rigorously established for n = 1 only, but we still may 
discuss the scaling limits at those values oi x. It is widely believed that at the 
critical values the model has a conformally invariant scaling limit. Moreover, the 
corresponding criticalities under renormalization are supposed to be unstable and 
stable correspondingly, so for x — x c there should be one conformally invariant 
scaling limit, whereas for the interval x € (x Cl oo) another, corresponding to x c . 
The scaling limit for low temperatures x S (0, x c ), a straight segment, is not 
conformally invariant. 

Plugging in n = 1 we obtain weight 

^length of loops 



Assigning the spins ±1 (represented by grey and white colors in Figure [T]) to sites 
of triangular lattice, we rewrite the weight as 

pairs of neighbors of opposite spins 

obtaining the Ising model (where the usual parameterization is exp(— 2/3) = x). 
The critical value is known to be /3 C = log 3/4, so one gets the Ising model at critical 
temperature for n = 1, x = l/\/3. A computer simulation of the Ising model on 
the square lattice at critical temperature, when the probability of configuration is 
proportional to ([1]) , is shown in Figure [2] 

For n = 1, x = 1 we obtain critical site percolation on triangular lattice. Taking 
n = (which amounts to considering configurations with no loops, just a curve 
running from a to b), one obtains for x c = l/V 2 + \/2 a version of the self-avoiding 
random walk. 

The following conjecture (see e.g. [T5]) is a direct consequence of physics pre- 
dictions and SLE calculations: 

Conjecture 3. For n £ [0,2] and x = x c (n), as lattice step goes to zero, the law 
of the interface converges to Schramm- Loewner Evolution with 



k = Am I (2ir — arccos(— n/2)) 
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Figure 2. The Ising model at critical temperature on a square. White and grey sites 
represent ±1 spins. Dobrushin boundary conditions (grey on lower and left sides, white 
on upper and right sides) produce, besides loop interfaces, an interface from the upper 
left to the lower right corner, pictured in black. When lattice step goes to zero, the law of 
the interface converges to SLE (3) , which is a conformally invariant random curve, almost 
surely simple and of Hausdorff dimension 11/8. 

For n £ [0,2] and x € (x c ,oo) (in particular for x = x c ), as lattice step goes to 
zero, the law of the interface converges to Schramm- Loewner Evolution with 

k = 4ir/ arccos(— n/2) . 

Note that to address this question one does not need to prove that the Nienhuis 
temperature is indeed critical (Conjecture [2]). 

We discussed loops on the hexagonal lattice, since it is a trivalent graph and 
so at most one interface can pass through a vertex. One can engage in similar 
considerations on the square lattice with special regard to a possibility of two 
interfaces passing through the same vertex, in which case they can be split into 
loops in two different ways (with different configuration weights). In the case of 
Ising (n = 1) this poses less of a problem, since number of loops is not important. 
For n = 1 and x = 1 we get percolation model with p = 1/2, but for a general 
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lattice this p need not be critical, so e.g. critical site percolation on the square 
lattice does not fit directly into this framework. 

2.3. Fortuin-Kasteleyn random cluster models. Another interest- 
ing class is Fortuin-Kasteleyn models, which are random cluster representations 
of g-state Potts model. The random cluster measure on a graph (a piece of the 
square lattice in our case) is a probability measure on edge configurations (each 
edge is declared either open or closed), such that the probability of a configuration 
is proportional to 



where clusters are maximal subgraphs connected by open edges. The two param- 
eters are edge- weight p € [0, 1] and cluster- weight q e (0, oo), with q £ [0, 4] being 
interesting in our framework (similarly to the previous model, q > 4 exhibits dif- 
ferent behavior). For a square lattice (or in general any planar graph) to every 
configuration one can prescribe a cluster configuration on the dual graph, such that 
every open edge is intersected by a dual closed edge and vice versa. See Figured] 
for a picture of two dual configurations with respective open edges. It turns out 
that the probability of a dual configuration becomes proportional to 



with the dual to p value = p*(p) satisfying p*/(l — p*) = q(l — p)/p. For 
p = p a d := yfq/(y/q+ 1) the dual value coincides with the original one: one gets 
Psd = {psd)* and so the model is self-dual. It is conjectured that this is also the 
critical value of p, which was only proved for q = 1 (percolation) , q = 2 (Ising) 
and q > 25.72. 

Again we introduce Dobrushin boundary conditions: wired on the counter- 
clockwise arc ab (meaning that all edges along the arc are open) and dual-wired on 
the counterclockwise arc ba (meaning that all dual edges along the arc are open, or 
equivalently all primal edges orthogonal to the arc are closed) - see Figure |3] Then 
there is a unique interface running from a to 6, which separates cluster containing 
the arc ab from the dual cluster containing the arc ba. 

We will work with the loop representation, which is similar to that in 12.21 
The cluster configurations can be represented as Hamiltonian (i.e. including all 
edges) non-intersecting (more precisely, there are no "transversal" intersections) 
loop configurations on the medial lattice. The latter is a square lattice which 
has edge centers of the original lattice as vertices. The loops represent interfaces 
between cluster and dual clusters and turn by ± 5 at every vertex - see Figure G3 
It is well-known that probability of a configuration is proportional to 



p 



# open edges 



(i-p) 



# closed edges # clusters 



# dual open edges 



\# dual closed edges # dual clusters 
P* ) H i 




which for the self-dual value p — p s d simplifies to 



# loops 



(2) 
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Figure 3. Loop representation of the random cluster model. The sites of the original 
lattice are colored in black, while the sites of the dual lattice are colored in white. Clusters 
and dual clusters and loops separating them are pictured. Under Dobrushin boundary 
conditions besides a number of loops there is an interface running from a to b, which is 
drawn in bold. Weight of the configuration is proportional to i^/q) 1o °p s 



Dobrushin boundary conditions amount to introducing two vertices with odd num- 
ber of edges: a source o and a sink 6, which enforces a curve running form a to b 
(besides loops) - see Figure for a typical configuration. 

Conjecture 4. For all q € [0,4], as the lattice step goes to zero, the law of the 
interface converges to Schramm- Loewner Evolution with k — Air/ arccos(— y/q/2). 

Conjecture was proved by Greg Lawler, Oded Schramm and Wendelin Werner 
[23] for the case of q = 0, when they showed that the perimeter curve of the 
uniform spanning tree converges to SLE (8). Note that with Dobrushin boundary 
conditions loop representation still makes sense for q = 0. In fact, the formula 
([21) means that we restrict ourselves to configurations with no loops, just a curve 
running from a to b (which then necessarily passes through all the edges), and all 
configurations are equally probable. 

Below we will outline our proof [39] |40] that for the Ising parameter q = 2 the 
interface converges to SLE (16/3), see Figured] It almost directly translates into a 
proof that the interface of the spin cluster for the Ising model on the square lattice 
at the critical temperature (which can be rewritten as the loop model in 12.21 for 
n = 1, only on the square lattice) converges to SLE (3), as shown on Figure [2] It 
seems likely that it will work in the n = 1 case for the loop model on hexagonal 
lattice described above, providing convergence to SLE (3) for x = x c and (a new 
proof of) convergence to SLE (6) for x — x c (and possibly for all x > x c ). 

Summing it up, Conjecture [3] was proved earlier for n = 1, O [37], 
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Figure 4. Interface in the random cluster Ising model at critical temperature with Do- 
brushin boundary conditions (loops not pictured). The law converges to SLE (16/3) when 
mesh goes to zero, so in the limit it has Hausdorff dimension 5/3 and touches itself almost 
surely. The random cluster is obtained by deleting some bonds from the spin cluster, so 
the interfaces are naturally different. Indeed, they converge to different SLE's and have 
different dimensions. However they are related: conjecturally, the outer boundary of the 
(non-simple) pictured curve and the (simple) spin interface in Figure [2] have the same 
limit after appropriate conditioning. 

whereas Conjecture [4] was established for q — 0, see [23]. We outline a technique, 
which seems to prove conformal invariance in two new cases, and provide new 
proofs for the only cases known before, making Conjecture @] solved for q = 
and q = 2, and Conjecture [3] for n = 1. The method also contributes to our 
understanding of universality phenomenon. 

Much of the method works for general values of n and q. The most interesting 
values of the parameters (where it does not yet work all the way) are n = 0, 
related to the self-avoiding random walk, and q = 1, equivalent to the critical bond 
percolation on the square lattice (in the latter case some progress was achieved by 
Vincent Beffara by a different method). Hopefully the lemma (essentially the 
discrete analyticity statement - see below) required to transfer our proof to other 
models will be worked out someday, leading to full resolution of these Conjectures. 



Towards conformal invariance of 2D lattice models 



11 



3. Schramm-Loewner Evolution 

3.1. Loewner Evolution. Loewner Evolution is a differential equation for 
a Riemann uniformization map for a domain with a growing slit. It was introduced 
by Charles Loewner in [25 in his work on Bieberbach's conjecture. 

In the original work, Loewner considered slits growing towards interior point. 
Though such radial evolution (along with other possible setups) is also important 
in the context of lattice models and fits equally well into our framework, we will 
restrict ourselves to the chordal case, when the slit is growing towards a point on 
the boundary. 

In both cases we choose a particular Riemann map by fixing its value and 
derivative at the target point. Chordal Loewner Evolution describes uniformization 
for the upper half-plane C+ with a slit growing from to oo (one deals with a 
general domain f2 with boundary points a, b by mapping it to C+ so that a i— ► 0, 
b h- > oo). 

Loewner only considered slits given by smooth simple curves, but more gener- 
ally one allows any set which grows continuously in conformal metric when viewed 
from oo. We will omit the precise definition of allowed slits (more extensive dis- 
cussion in this context can be found in [22]), only noting that all simple curves are 
included. The random curves arising from lattice models (e.g. cluster perimeters 
or interfaces) are simple (or can be made simple by altering them on the local 
scale) . Their scaling limits are not necessarily simple, but they have no "transver- 
sal" self-intersections. For such a curve to be an allowed slit it is sufficient if it 
touches itself to never venture into the created loop. This property would follow 
if e.g. a curve visits no point thrice. 

Parameterizing the slit 7 in some way by time t, we denote by gt{z) the con- 
formal map sending C+ \ 7* (or rather its component at 00) to C+ normalized 
so that at infinity gt(z) = z + a(t)/z + C(l/|z| 2 ), the so called hydrodynamic 
normalization. It turns out that a(t) is a continuous strictly increasing function 
(it is a sort of capacity- type parameter for 74), so one can change the time so that 



Denote by w(t) the image of the tip j(t). The family of maps gt (also called a 
Loewner chain) is uniquely determined by the real-valued "driving term" w(t). 
The general Loewner theorem can be roughly stated as follows: 

Loewner's theorem. There is a bisection between allowed slits and continuous 
real valued junctions w(t) given by the ordinary differential equation 




(3) 



dtgt(z) 



2 



9o{z) 



(4) 



9t(z) 



w(t) 



= z 



The original Loewner equation is different since he worked with smooth radial 
slits and evolved them in another (but related) way. 



12 



Stanislav Smirnov 



3.2. Schramm-Loewner Evolution. While a deterministic curve 7 cor- 
responds to a deterministic driving term w(t), a random 7 corresponds to a random 
w(t). One obtains SLE (k) by taking w(t) to be a Brownian motion with speed k: 

Definition 5. Schramm-Loewner Evolution, or SLE (k), is the Loewner chain one 
obtains by taking w(t) = ^/KB t , K € [0, 00). Here B t denotes the standard (speed 
one) Brownian motion (Wiener process). 

The resulting slit will be almost surely a continuous curve. So we will also 
use the term SLE for the resulting random curve, i.e. a probability measure on 
the space of curves (to be rigorous one can think of a Borel measure on the space 
of curves with uniform norm). Different speeds k produce different curves: we 
grow the slit with constant speed (measured by capacity), while the driving term 
"wiggles" faster. Naturally, the curves become more "fractal" as k increases: for 
k < 4 the curve is almost surely simple, for 4 < n < 8 it almost surely touches 
itself, and for n > 8 it is almost surely space-filling (i.e. visits every point in 
C+) - see [55J 130] for these and other properties. Moreover, Vincent Beffara [7] 
has proved that the Hausdorff dimension of the SLE (k) curve is almost surely 
min(l + k/8,2). 

3.3. Confer mal Markov property. Suppose we want to describe the 
scaling limits of cluster perimeters, or interfaces for lattice models assuming their 
existence and conformal invariance. We follow Oded Schramm [35] to show that 
Brownian motion as the driving force arises naturally. Consider a simply connected 
domain O with two boundary points, a and b. Superimpose a lattice with mesh 
e and consider some lattice model, say critical percolation with the Dobrushin 
boundary conditions, leading to an interface running from a to b, which is illus- 
trated by Fig. [1] for a rectangle with two opposite corners as a and b. So we end up 
with a random simple curve (a broken line) connecting a to b inside Q. The law of 
the curve depends of course on the lattice superimposed. If we believe the physi- 
cists' predictions, as mesh tends to zero, this measure on broken lines converges 
(in an appropriate weak-* topology) to some measure /1 = /i(fi, a, b) on continuous 
curves from a to b inside Q. 

In this setup the conformal invariance prediction can be formulated as follows: 

(A) Conformal Invariance. For a conformal map 4> of the domain f2 one has 

(/*(«, a, 6)) = M (0(0), <K<*), <K&)) • 

Here a bijective map : f2 — > (f>(Q) induces a map acting on the curves in Q, which 
in turn induces a map on the probability measures on the space of such curves, 
which we denote by the same letter. By a conformal map we understand a bijection 
which locally preserves angles. 

Moreover, if we start drawing the interface from the point a, we will be walking 
around the grey cluster following the right-hand rule - see Figure 1. If we stop at 
some point a' after drawing the part 7' of the interface, we cannot distinguish the 
boundary of Q from the part of the interface we have drawn: they both are colored 
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Figure 5. (A) Conformal Invariance: conformal image of the law of the curve 7 
(dotted) in Q coincides with the law of the curve 7 in the image domain (j>(Sl). 



grey on the (counterclockwise) arc a'b and white on the arc ba' of the domain 
So we can say that the conditional law of the interface (conditioned on it starting 
as 7') is the same as the law in a new domain with a slit. We expect the limit law 
\x to have the same property: 

(B) Markov property. The law conditioned on the interface already drawn is 
the same as the law in the slit domain: 

n(n,a,b)W = n(nw,a',b) . 



identical laws 




Figure 6. (B) Markov property: The law conditioned on the curve already drawn is 
the same as the law in the slit domain. In other words when drawing the curve we do 
not distinguish its past from the boundary. 

If one wants to utilize these properties to characterize /i, by (A) it is sufficient 
to study some reference domain (to which all others can be conformally mapped), 
say the upper half-plane C+ with a curve running from to 00. Given (A), the 
second property (B) is easily seen to be equivalent to the following: 
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(B') Conformal Markov property. The law conditioned on the interface al- 
ready drawn is a conformal image of the original law. Namely, for any conformal 
map G — Gy from C + \ 7' to C + preserving 00 and sending the tip of 7' to 0, we 
have 

/x(C + ,0,oo)| 7 ' - G- 1 ( M (C+,0,oo)) . 



identical laws 

Q-l / \ 

m(C+,0,co) .— > G- 1 (MC+,0,oo)) MC + ,0,oo)|7' 




Figure 7. (B') Conformal Markov property: The law conditioned on the curve 
already drawn is a conformal image of the original law. In other words the curve has 
"Markov property in conformal coordinates" . 

Remark 6. Note that property (B') is formulated for the law fj, for one domain 
only, say C + as above. If we extend fi to other domains by conformal maps, it 
turns out that (A) and (B) are equivalent to (B') together with scale invariance 
(under maps z 1— > kz, k > 0). 

To use the property (B'), we describe the random curve by the Loewner Evolu- 
tion with a certain random driving force w(t) (we assume that the curve is almost 
surely an allowed slit). If we fix the time t, the property (B') with the slit 7[0,i] 
and the map G t {z) = g t {z) — w(t) can be rewritten for random conformal map 
G t +8 conditioned on G t (which is the same as conditioning on 7[0,i]) as 

G t +s\G t = G t (Gs) ■ 
Expanding G's near infinity we obtain 

z-w(t + S) + . ..\G t = (z- w(t) + ... )o(z - w(S) + ...) = z-(w(t)+w(5)) + . . . , 
concluding that 

w(t + S) - w(t)\G t = w(S) . 

This means that w(t) is a continuous (by Loewner's theorem) stochastic process 
with independent stationary increments. Thus for a random curve satisfying (B') 
the driving force w(t) has to be a Brownian motion with a certain speed k e [0, 00) 
and drift ael: 

w(t) — \fnB t + at . 
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Applying (A) with anti-conformal reflection <f>(u + iv) = —u + iv or with stretching 
<fi(z) — 2z shows that a vanishes. So one logically arrives at the definition of 
SLE and the following 

Schramm's principle. A random curve satisfies (A) and (B) if and only if it is 
given by SLE (n) for some n € [0,oo). 

The discussion above is essentially contained in Oded Schramm's paper [32] for 
the radial version, when slit is growing towards a point inside and the Loewner dif- 
ferential equation takes a slightly different form. To make this principle a rigorous 
statement, one has to require the curve to be almost surely an allowed slit. 



4. SLE as a scaling limit 

4.1. Strategy. In order to use the above principle one still has to show the 
existence and conformal invariance of the scaling limit, and then calculate some 
observable to pin down the value of n. For percolation one can employ its locality 
or Cardy's formula for crossing probabilities to show that k — 6. Based on this 
observation Oded Schramm concluded in |32] that if percolation interface has a 
conformally invariant scaling limit, it must be SLE (6). 

But it is probably difficult to show that some interface has a conformally in- 
variant scaling limit without actually identifying the latter. 

To identify a random curve in principle one needs "infinitely many observables," 
e.g. knowing for any finite number of points the probability of passing above them. 
This seems to be a difficult task, which is doable for percolation since the locality 
allows us to create many observables from just one (crossing probability), see [38) . 

Fortunately it turns out that even in the general case if an observable has a 
limit satisfying analogues of (A) and (B), one can deduce convergence to SLE (k) 
(with k determined by the values of the observable) . 

This was demonstrated by Greg Lawler, Oded Schramm and Wendelin Werner 
in [23] in establishing the convergence of two related models: of loop erased random 
walk to SLE (2) and of uniform spanning tree to SLE (8). 

They described the discrete curve by a Loewner Evolution with unknown ran- 
dom driving force. Stopping the evolution at times t and s and comparing the 
values of the observable, one deduces (approximate) formulae for the conditional 
expectation and variance of the increments of the driving force. Skorokhod embed- 
ding theorem is then used to show that driving force converges to the Brownian 
motion. Finally one has to prove a (stronger) convergence of the measures on 
curves. The trick is that knowing just one observable (but for all domains) after 
conditioning translates to a continuum of information about the driving force. 

We describe a different approach with the same general idea, which is perhaps 
more transparent, separating "exact calculations" from "a priori estimates". The 
idea is to first get a priori estimates, which imply that collection of laws is precom- 
pact in a suitable space of allowed slits. Then to establish the convergence it is 
enough to show that limit of any converging subsequence coincides with SLE. To 
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do that wc describe the subscqucntial limit by Loewner Evolution (with unknown 
random driving force w(t)) and extract from the observable enough information to 
evaluate expectation and quadratic variation of increments of w(t). Levy's charac- 
terization implies that w(t) is the Brownian motion with a particular speed k and 
so our curves converge to SLE (k). 

As an example we discuss below an alternative proof of convergence to SLE (6) 
in the case of percolation, which uses crossing probability as an observable. For 
a domain Q with boundary points a, b superimpose triangular lattice with mesh 
e and Dobrushin boundary conditions. Wc obtain an interface j e (between open 
vertices on one side and closed on another) running from a to b, see Fig. [1] i.e. a 
measure /i e on random curves running from a to b. 

4.2. Compactness. First we note that the collection {fi e } is precompact (in 
weak-* topology) in the space of continuous curves that are Loewner allowed slits. 

The necessary framework for precompactness in the space of continuous curves 
was suggested by Michael Aizenman and Almut Burchard [5] . It turns out that ap- 
propriate bounds for probability of an annulus being traversed k times imply tight- 
ness: a curve has a Holder parameterization with stochastically bounded norm. 
Hence {fJ. e } is precompact by Prokhorov's theorem: a (uniformly controlled) part 
of [i t is supported on a compact set (of curves with norm bounded by M), and so 
such parts are weakly precompact by Banach-Alaoglu theorem, whereas the mass 
of the remainder tends uniformly to zero as M — > oo. 

The curves on the lattice are simple, so they cannot have transversal self- 
intersections even after passing to the limit. So to check that for any weak limit of 
/U e 's almost every curve is an allowed slit, one has to check that as we grow it the 
tip is always visible and moves continuously when viewed from infinity. Essentially, 
one has to rule out two scenarios: that the curve passes for a while inside already 
visited set, and that the curve closes a loop, and then travels inside before exiting. 
Both are reduced to probabilities of annuli traversing. 

In the case of percolation one uses the Russo-Seymour- Welsh theory [3lJ [36] 
together with Michael Aizenman's observation pQ (that in the limit interface can 
visit no point thrice - "no 6 arms") to obtain the required estimates. 

Since the collection of interface laws {/i £ } is precompact (in weak-* topology) 
in the space of continuous curves that are Loewner allowed slits, to show that as 
mesh goes to zero the interface law converge to the law of SLE (6), it is sufficient 
to show that the limit of any converging subsequence is in fact SLE (6). 

Take some subsequence converging to a random curve in the domain fl from a 
to b. We map conformally to a half-plane C+, obtaining a curve 7 from to 00 
with law fi. We must show that fi is given by SLE (6). 

By a priori estimates 7 is almost surely an allowed slit. So we can describe 7 
by a Loewner Evolution with a (random) driving force w(t). It remains to show 
that w(t) — y6Bt. Note that at this point we only know that w(t) is an almost 
surely continuous random function - we do not even have a Markov property. 
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4.3. Martingale observable. Given a topological rectangle (a simply con- 
nected domain Q with boundary points a, b, c, d) one can superimpose a lattice 
with mesh e onto fl and study the probability II £ (ft, [a,b], [c,d]) that there is an 
open cluster joining the arc [a, b] to the arc [c, d] on the boundary of ft. It is 
conjectured that there is a limit II := limg^oITe, which is conformally invariant 
(depends only on the conformal modulus of the configuration O, a, b, c, d), and 
satisfies Cardy's formula (predicted by John Cardy in [9] and proved in [37) ) in 
half-plane: 

ncc,, [!-«,!], [ceo]) = r(1 ^ /3) » 1/3 ^Q.|;|;») == *■(«)■ (s) 

Above 2-F1 is the hypergeometric function, so one can alternatively write 

F(u) = (v(l - v))~ 2/3 dv I J (v{l -v))~ 2/3 dv . 

Particular nature of the function is not important, we rather use the fact that there 
is an explicit formula for half-plane with four marked boundary points and hence 
by conformal invariance for an arbitrary topological rectangle. The value k = 6 
will arise later from some expression involving derivatives of F. 

Assume that for some percolation model we are able to prove the above conjec- 
ture (for critical site percolation on the triangular lattice it was proved in [3"Tll38| ). 

Add two points on the boundary, making f2 a topological rectangle axby and 
consider the crossing probability II £ (fi, [a, x], [b, y]) (from the arc ax to the arc by 
on a lattice with mesh e). 

Parameterize the interface j e in some way by time, and draw the part 7 £ [0,t]. 
Note that it has open vertices on one side (arc y e (t)a) and closed on another (arc 
wf e (t)). Then any open crossing from the arc by to the arc ax inside ft is either 
disjoint from 7 £ [0,i], or hits its "open" arc "f e (t)a. In either case it produces an 
open crossing from the arc by to the arc J e (t)x inside fl \ j e [0,t], and converse 
also holds. Therefore one sees that for every realization of y e [0, t] the crossing 
probability conditioned on 7 e [0,£] coincides with crossing probability in the slit 
domain ft \ 7 e [0, t]: 

n e (ft, [a,x], [b, y]| 7e [0, t]) = n £ (fl \ 7e [0, t], [ 7e (t), x], [b, y]) , (6) 

- an analogue of the Markov property (B). Alternatively this follows from the fact 
that II can be understood in terms of the interface as the probability that it touches 
the arc xb before the arc by. For example, on Figure [T] there is no horizontal grey 
(open) crossing (there is a vertical white crossing instead), and interface traced 
from the left upper corner a touches the lower side xb before the right side by. 

Stopping the curve at times t < s and using ^ we can write by the total 
probability theorem for every realization of 7 e [0,t] 

n £ (n\7 e [o,t],[7 e (t), a], [&,»]) 

= E 7eM (n £ (ft\ 7£ [0, S ],[ 7£ ( S ), a ;],[6, 2 /])|7[0,t]) . (7) 
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The same a priori estimates as in the previous subsection show that the identity 
||7J) also holds for the (subsequential) scaling limit \i (strictly speaking there is an 
error term in case the interface touches the arcs [ax] or [ya] before time s, but it 
decays very fast as we move x and y away from a). We know that the scaling limit 
II := lim e ^o n e of the crossing probabilities exists and is conformally invariant, so 
we can rewrite ([7]) for the curve 7 with Loewner parameterization as 



n(C + \7[0,t],[7(t),z],[cx> ) tf]) 

= E 7[M] (II (C+ \ 7 [0, s], h(s),x], [00, y]) | 7 [0, t\) 



(8) 



for almost every realization of 7 [0,i]. Moreover we can plug in exact values of 
the crossing probabilities, given by the Cardy's formula. Recall that the domain 
C+ \ 7 [0,t] is mapped to half-plane by the map gt{z) with 7 (t) 1— * w(t). Then the 
map z i-» llfylfX) also maps it to half-plane with j(t) i-> g^-gtiv) ' y ^ °' 
DC 1 — >• 1. Using conformal invariance and applying Cardy's formula we write 



n(c + \ 7 [o,t],[ 7 W,x],[oo, 2/ ]) = n(c 

= F 



9t{y) - w(t) 
9t(x)-g t (y)' 

g t (x) - w(t) 



, [oo,0] 



9t(x)-g t (y)J 



(9) 



for Cardy's hypergeometric function F. 



4.4. Conformally invariant martingale. Plugging © into both sides 



of 



we arrive at 



F 



g t {x) - w(t) 
9t(x) - g t (y) 



= E 



7[M] 



F 



g s (x) - w{s) 

g s (x)-9s(y) 



h[o,t] 



(10) 



Remark 7. Denote by x t := gt(%) — w(t) and y t :— gt(y) — w(t) trajectories 
of x and y under the random Loewner flow. Then (j 1 Of) essentially means that 

F y x ^l yt ^ is a martingale. 



Since we want to extract the information about w(t), we fix the ratio x/{x — 
y) :— 1/3 (anything not equal to 1/2 would do) and let x tend to infinity: y :— —2x, 
x —y +00. Using the normalization gt(z) = z + 2t/z + 0(l/z 2 ) at infinity, writing 
Taylor expansion for F, and plugging in values of derivatives of F at 1/3, we obtain 
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the following expansion for the right-hand side of pop : 

x- w(t) + 2t/x + 0(l/x 2 ) 



= F 



(x + 2t/x + 0(l/x 2 )) - (-2x + 2t/(-2x) + 0{l/x 2 )) 



3 3 a; 3 x 2 \x 3 

r 1 1) - (I) I + (If- (I) + gtr iTlU + o ' 1 



3 3 V 3 / a; V 3 V 3 J 3 2 • 2 V 3 



r(2/3) 



^K(^-S)+o(l 



x 2 r(l/3)r(4/3) 3 2 / 3 
=: A - -BEw(t) - \cE (w(t) 2 -6t)+0 ( 

where we plugged in values of the derivative for hypergeometric function. Using 
similar reasoning for the right-hand side of (|10p we arrive at the following identity: 



—BEw{t) - —,C E (w{i) 2 -6t)+0 ( ^ 

X X \ X 



= A- ^BE 7[fiS] (w(sM0,t]) - ^CE l[Ls] (w(s) 2 - 6s| 7 [0,i]) +0 

Equating coefficients in the series above, we conclude that 

E w[t , s] (w(s)\w[0, t])=0, E w[t<s] (w{s) 2 - 6s\w[0, t}) = w{tf - 6t . (11) 

Thus w(t) is a continuous (by Loewner's theorem) process such that both 

w{t) and w(t) 2 - 6t 

are martingales so by Levy's characterization of the Brownian motion w(t) = 
V6B t , and therefore SLE (6) is the scaling limit of the critical percolation interface. 

The argument will work wherever Cardy's formula and a priori estimates are 
available, particularly for triangular lattice. More generally, any conformally in- 
variant martingale will do, with value of k arising from its Taylor expansion. 

Remark 8. The scheme can also be reversed to do calculations for SLE's, if an 
observable is a martingale (e.g. crossing probability). Indeed, writing the same 
formulae with x/(x — y) = a we conclude that the coefficient by 1/x 2 , namely 

C = M y 2a) tF'(a) + ^E (w(t) 2 )F"(a) , 
1 — a 2 
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vanishes. Since for w(t) = \^6B(t) one has E(u>(i) 2 ) = 6t, we arrive at the 
differential equation 

f-%'( a)+ F>) = C, 
6(1 — a) 

With the given boundary data it has a unique solution, which is Cardy's hyperge- 
ometric function. 



5. Ising model and beyond 

The martingale method as described above shows that to construct a conformally 
invariant scaling limit for some model we need a priori estimates and a non-trivial 
martingale observable with a conformally invariant scaling limit. 

5.1. A priori estimates. A priori estimates are necessary to show that 
collection of interface laws is precompact in weak-* topology (on the space of 
measures on continuous curves which are allowed slits). 

If we follow the same route as for percolation (via the work [2] of Michael Aizen- 
man and Almut Burchard), we only need to evaluate probabilities of traversals of 
an annulus in terms of its modulus. For percolation such estimates are (almost) 
readily available from the Russo-Seymour- Welsh theory. For uniform spanning tree 
and loop erased random walk one can derive the estimates using random walk con- 
nection and the known estimates for the latter (a "branch" of a uniform spanning 
tree is a loop erased random walk) , see [3J [35] . 

For the Ising model the required estimates do not seem to be readily available, 
but a vast arsenal of methods is at hand. Essentially all we need can be reduced 
by monotonicity arguments to spin correlation estimates of Bruria Kaufman, Lars 
Onsager and Chen Ning Yang [HI Wf\ . 

For general random cluster or loop models such exact results are not available, 
but we actually need much weaker statements, and many of the techniques used 
by us for the Ising model (like FKG inequalities) arc well-known in the general 
case. 

So this part does not seem to be the main obstacle to construction of scaling 
limits, though it might require very hard work. Moreover, following the proposed 
approach we actually get that interfaces have a Holder parameterization with uni- 
formly stochastically bounded norm. Thus rather weak kinds of convergence of 
interfaces would lead to convergence in uniform norm (or rather weak-* conver- 
gence of measures on curves with uniform norm) . 

It also appears that the same a priori estimates can be employed to show 
observable convergence in the cases concerned, and hopefully they will be sufficient 
for other models. So a more pressing question is how to construct a martingale 
observable. 

5.2. Conformally covariant martingales. Suppose that for every sim- 
ply connected domain Q with a boundary point a we have defined a random curve 
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7 starting from a. Mark several points 6, c, . . . in f2 or on the boundary. Remark [7J 
suggests the following definition: 

Definition 9. We say that a function (or rather a differential) F(Sl, a,b,c, . . . ) is 
a conformal (covariant) martingale for a random curve 7 if 

F is conformally covariant : F (0, a, b, c, . . . ) = 
=F (0(0), 0(a), 0(6), 0(c), . . . ) • 0'(&)V(&) V(c)V(c) 5 . . . , (12) 

and 

F(fi\7[0,t],7(t),&,c,...) is a martingale (13) 
with respect to the random curve 7 drawn from a (with Loewner parameterization). 

Introducing covariance at 6, c, . . . we do not ask for covariance at a, since it 
always can be rewritten as covariance at other points. And applying factor at 
a would be troublesome: once we started drawing a curve the domain becomes 
non-smooth in its neighborhood, creating problems with the definition. 

If the exponents a, (3, . . . vanish, we obtain an invariant quantity. While the 
crossing probability for the percolation was invariant, many quantities of interest in 
physics are covariant differentials, e.g. cluster density at c would scale as a lattice 
step to some power (depending on the model), so we would arrive at a factor 

|0'( C )| 5 =0'( C ) 5 / 2 0'( C ) A "/2 . 

There are other possible generalizations, e.g. one can add the Schwarzian derivative 
of to (I12|l . The two properties in Definition [9] are analogues of (A) and (B), and 
similarly combined they show that for the curve 7 mapped to half-plane from any 
domain f2 so that a 1— ► 0, b > 00, c 1— > x (note that the image curve in C+ might 
depend on f2 - we only know the conformal invariance of an observable, not of 
the curve itself) we have an analogue of (B'), which was already mentioned in 
Remark [7] for percolation. Namely 

F (C + ,0,^,g t (x), . . .) ■ g'(xrg' t (x) 5 . . . , 

is a martingale with respect to the random Loewner evolution (covariance factor 
at b = 00 is absent, since g' t (oo) = 1). 

The equation pop can be written for this F, and if we can evaluate F exactly, 
the same machinery as one used by Greg Lawler, Oded Schramm and Wendelin 
Werner in |23j or as the one discussed above for percolation proves that our random 
curve is SLE. So one arrives at a following generalization of Oded Schramm's 
principle: 

Martingale principle. If a random curve 7 admits a (non-trivial) conformal 
martingale F, then 7 is given by SLE with k (and drift depending on modulus of 
the configuration) derived from F. 
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Remark 10. In chordal situation we consider curves growing from a towards 
another boundary point & in a simply connected domain. But the same conclusion 
would hold on general domains or Riemann surfaces with boundary once we hnd 
a covariant martingale (for appropriate generalizations of Loewner Evolutions see 
e.g. the book [4]). The only difference is that driving force of the corresponding 
Loewner Evolution will be a Brownian motion with drift depending on conformal 
modulus of the configuration Q, a, b, c, . . . , leading to SLE generalizations. Starting 
from lattice models with various boundary conditions and conditioned on various 
events, one can see which drifts will be of interest for SLE generalizations. 

5.3. Discrete analyticity. Passing to the lattice model, we want to find a 
discrete object, which in the limit becomes a conformally covariant martingale. 

Martingale property is actually more accessible in the discrete setting. For 
example, functions which are defined as observables (like probability of the interface 
going through a vertex, edge density for the model, etc.) have the martingale 
property built in, and so only conformal covariance must be established. 

Alternatively, one can work with a discrete function _F(f2, a, 6, c) (a priori not 
related to lattice models) which has a conformally covariant scaling limit by con- 
struction. Then we need to connect it to a particular lattice model, establishing a 
martingale property (fTU|) . In the discrete case it is sufficient to check the latter for 
a curve advanced by one step. Assume that once we have drawn the part 7' of the 
interface from the point a to point a', it turns left with probability p = p(f2, 7' ,a! ', b) 
creating a curve 7; = 7 U {a{\ or right with probability (1 — p) creating a curve 
7 r = 7 U {a r }. Then it is enough to check the identity 

F(J1 \ 7', a', b, c) = P F(n \ lh a h b, c) + (1 - p)F(Sl \ 7r , a r , b, c) , (14) 

p = p(Sl, 7', a', b) . 

Actually our proof for the Ising model can be rewritten that way, with F defined 
as a solution of an appropriate discretization of the Riemann Boundary Value 
Problem (jTTJ) - the observable nature of F never comes up. 

Moreover, starting with F one can define a random curve by choosing "turning 
probabilities" p so that identity (I14|) is satisfied, obtaining a model with confor- 
mally invariant scaling limit by "reverse engineering." For example, starting with 
a harmonic function of c with boundary values 1 on the arc ba and on the arc 
ab, one obtains a unique discrete random curve, which has it as a martingale. 
Note that such a function is a particular case a = 1 of the martingale (fT5|) below, 
corresponding to k = 4 (or rather its integral). In [35] Oded Schramm and Scott 
Sheffield introduced this curve with a nicer "Harmonic Explorer" definition, and 
utilizing the mentioned observable showed that it indeed converges to SLE (4) . It 
seems that in this way one can use the solutions to the problem (jTTJ) to construct 
models converging to arbitrary SLE's, however it is not clear though whether they 
would similarly have "nicer" definitions. 

Anyway, for either approach to work we need a discrete conformal covariant 
with a scaling limit. We have tried discretizations of many conformally invariant 
objects (extremal length, capacity, solutions to variational problems, . . . ) and the 
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most promising in this context seem to be discrete harmonic or analytic functions 
in additional variable(s) (in c, . . . ). Firstly, all other invariants can be rewritten in 
this way. Secondly, discretization of harmonic and analytic functions is a nice and 
very well studied (especially in the case of harmonic ones) object. Thirdly, one can 
obtain very non-trivial invariants by just checking local conditions: harmonicity or 
analyticity inside plus some boundary conditions (Dirichlet, Neumann, Riemann- 
Hilbert, etc.). The most natural candidate would be a harmonic function solving 
some Dirichlet problem. 

Note that such an observable is known for the Brownian motion. A classical 
theorem [13] of Shizuo Kakutani states that in a domain f2 exit probabilities for 
Brownian motion started at z are harmonic functions in z with easily determinable 
boundary values. Though Kakutani works directly with Brownian motion, one can 
do the same for the random walk (which is actually much easier, since discrete 
Laplacian of the exit probability is trivially zero), and then passing to a limit 
deduce statements about Brownian motion, including its conformal invariance. 

5.4. Classification of conformal martingales. Before we start work- 
ing in the discrete setting, we might want to investigate which functions are con- 
formal martingales for SLE curves, and so can arise as scaling limits of martingale 
observables for lattice models. 

As discussed in Remark [8j one can write partial differential equations for 
SLE conformal martingales. For small number of points those equations can be 
solved, and in such a way one computes dimensions, scaling exponents and other 
quantities of interest. For any particular value of k we can see which martingales 
have the simplest form and so are probably easier to work with. Also if they have 
a geometric SLE interpretation (like probability of SLE curve going to one side of 
a point, etc.) we can study similar quantities for the lattice model. 

It turns out that only for k — 4 one obtains a nice harmonic martingale with 
Dirichlet boundary conditions. In that case the probability of SLE (4) passing to 
one side of a point z is harmonic in z and has boundary conditions and 1, see 
[33] , Oded Schramm and Scott Sheffield [35] constructed a model which has this 
property on discrete level built in. Unfortunately the property was not yet observed 
in any of the classical models conjecturally converging to SLE (4), though results 
of Kenyon |16] show it holds for double-domino curves in Temperley domains (i.e. 
a domain with the boundary satisfying a certain local condition). 

In the case of mixed Dirichlet-Neumann conditions, it becomes possible to 
work with some other values of k, including uniform spanning tree k = 8, which is 
exploited in [23]. There are also covariant candidates for a few other values of n 
(notably 8/3 which corresponds to self-avoiding random walk), but they were not 
yet observed in lattice models. 

Thus to study general models, one is forced to utilize more general boundary 
value problems with a Riemann(-Hilbert) Boundary Value Problem being the nat- 
ural candidate. Besides harmonic function it involves its harmonic conjugate, and 
so is better formulated in terms of analytic functions. Moreover, discrete analyt- 
icity involves a first order Cauchy-Riemann operator, rather than a second order 
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Laplacian, and so it should be easier to deal with than harmonicity. 

As discussed above, we can classify all analytic martingales. For chordal 
SLE and F with three points a, 6, z as parameters we discover two particularly nice 
families. The following Proposition will be discussed in [40l[42] and our subsequent 
work: 

Proposition 11. Let CI be a simply connected domain with boundary points a, b. 
Let &(z) — a, b, z) be a mapping offl to a horizontal strip R x [0, 1], such that 
a and b are mapped to ^foo. Then 

F(n, a, b, z) = $'(»" with a = - - 1 , (15) 

K 

is a martingale for SLE(k). Let ^f(z) — x f'(f2, a, 6, z) be a mapping of to a 
half-plane C+, such that a and b are mapped to oo and correspondingly. Then 

G{Q, a, b, z) = t>'(z) a y'{b)- a with «=---, (16) 

is a martingale for SLE (k) . 

These martingales make most sense for k € [4,8] and k G [8/3,8] correspond- 
ingly, and are related to observables of interest in Conformal Field Theory (which 
was part of our motivation to introduce them) . Note that both functions are co- 
variant with power a (which is the spin in physics terminology), and solve the 
Riemann boundary value problems 

Im(F(z)T(z) Q ) = 0, zedtt, (17) 

where t(z) is the tangent vector to dQ at z. 

The problem is to observe these functions in the discrete setting, and some 
intuition can be obtained from their geometric meaning for SLE's. For example, F 
is roughly speaking (one has to consider an intermediate scale to make sense of it) 
an expectation of SLE curve passing through z taken with some complex weight 
depending on the winding. 

5.5. Height models and Coulomb gas. The above-mentioned expec- 
tation actually makes more sense (and is immediately well-defined) in the discrete 
setting and one arrives at the same object with the same complex weight via several 
different approaches. 

One way is to consider the Coulomb gas arguments (cf. by Bernard Nien- 
huis) for the loop representation. In the random cluster case at criticality, the 
weight of a loop is ^fq - recall © . We randomly and independently orient the loops, 
and introduce the height function h which whenever a loop is crossed changes by 
±1 (depending on loop direction - think of a topographic map). One could weight 
oriented loops by y/q/2, obtaining essentially the same model. However it makes 
sense to consider a complex weight instead. When q is in the [0,4] range, there is 
a complex unit number fi = exp(fc • 2iri) such that 

k = — arccos (y/q/2) or fi + fl=^/q. (18) 
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We (independently and randomly) orient all loops, prescribing weight fj, per coun- 
terclockwise and /2 per clockwise loop. 

Forgetting orientation of loops reconstructs the original model. Unfortunately 
the new partition function is complex and no longer leads to a probability measure 
(moreover, its variation blows up as the lattice step goes to zero), but it can be 
defined locally, making it much more accessible. 

Indeed, going around a cycle, and turning by A z at vertex z, the total sum of 
turns X^ecycic ^ ^ s ±27r depending on whether the cycle is counter or clockwise. 
So the weight per cycle can be written as Jlzecycic ex P(*^ • A z ) and the total weight 
of the configuration is II z esi exp(ifc • A z ), which can be computed locally (without 
reference to the global order of cycles). The same weight can also be written in 
terms of the gradient of height function. 

The interface is always oriented from a to 6, so that the height function is always 
equal to on the arc ab and to 1 on the arc ba. From physics arguments the inter- 
face curve (being "attached" to the boundary on both sides) should be weighted 
differently from loops, namely by exp(i(2/c — 1/2) • A z ) per turn. When interface 
runs between two boundary points (being oriented from a to 6), these factors do 
not matter, since total turn from a to b is independent of the configuration. 

However, if we choose a point z on an interface and reverse the orientation of 
one of its halves (so that it is oriented from a to z and from z to &), the interface 
inputs a non-constant complex factor. This orientation reversal has a nice meaning: 
after it the height function acquires a +2 monodromy at z: when we go around z 
we cross two curves (halves of the interface) incoming into it. 

All the loops (when we forget their orientation) still contribute the same yfq per 
loop, and the complex weight can be expressed in terms of the interface winding 
(total turn expressed in radians) from b to z, denoted by w(j,b — > z). So one 
logically arrives at the partition function Z for our model with +2 monodromy at 
z: 

F(fi, a, 6, z) := Z +2 monodromy at z = ^Xze~< exp (i(4fc - l)w(j, b-> z)). (19) 

This function is clearly a martingale, and there are strong indications (both from 
mathematics and physics points of view) that it is discrete analytic. 

This follows from the fact that the interface can arrive at a boundary point z 
from b with a unique winding equal to the winding of the boundary from b to z, 
so we can express it in terms of the tangent vector r(z). Writing this down, we 
discover that the function F solves a discrete version of the Riemann Boundary 
Value Problem (JTTj) with a = 1 - 4fc. 

Remark 12. The continuum problem was solved by the function (fT5")) . so if we 
establish discrete analyticity it only remains to show that a solution to a discrete 
Riemann Boundary Value Problem converges to its continuum counterpart. More- 
over, combining identity a = 1 — 4fc with (fT5)) and (fTS)) we obtain the relation 
between k and q stated in Conjecture 2J 

This convergence problem seems to be difficult (and open in the general case). 
The way we solve it in the Ising case is sketched below. 
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Figure 8. Rearrangement at a point z: we only change connections inside a small circle 
marking z. Either interface does not pass through z in both configurations, or it passes 
in a way similar to the pictured above. On the left the interface (in bold) passes through 
z twice, on the right (after the rearrangement) it passes once, but a new loop through z 
appears (also in bold). The loops not passing through z remain the same, so the weights 
of configurations differ by a factor of yfq because of the additional loop on the right. 
To get some linear relation on values of F, it is enough to check that any pair of such 
configurations makes equal contributions to two sides of the relation. 



There are other indications that this function is nice to work with. Indeed, the 
easiest form of discrete analyticity involves local partial difference relations, and 
to prove those we should count configurations included into our expectation. To 
obtain relations, we need some bijections in the configuration space, and the easiest 
ones are given by local rearrangements (we worked with global rearrangements for 
percolation [37], but such work must be more difficult for non-local models). 

The easiest rearrangement involves redirecting curves passing through z, see 
Figure [SI and we have a good control over relative weights of configurations when- 
ever they are defined through windings. Counting how much a pair of configura- 
tions contributes to values of F at neighbors of z, we get some relations. Moreover, 
a careful analysis shows that the maximal number of relations is attained with the 
complex weight ([T9|) . 



5.6. Ising model. Wc finish with a sketch of our proof for the random cluster 
representation of the Ising model (i.e. q = 2) on the square lattice eZ 2 at the critical 
temperature. As before consider loop representation in a simply connected domain 
with two boundary points a and b and Dobrushin boundary conditions. 

Consider function F = F e (Cl,a,b, z) given by (| 19|) which is the expectation 
that interface from a to b passes through a vertex z taken with appropriate unit 
complex weight. Note that for Ising q — 2, so k — 1/8 and the weight is Fermionic 
(which of course was expected): a passage in the same direction but with a 27r 
twist has a relative weight —1, whereas a passage in the opposite direction with a 
counterclockwise 7r twist has a relative weight — i. 

As discussed F automatically has the martingale property when we draw 7 
starting from a, so only conformal invariance in the limit has to be checked. 
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Color lattice vertices in chessboard fashion, and to each edge e prescribe ori- 
entation such that it points from a black vertex to a white one, turning it into 
a vector, or equivalently a complex number e. Denote by 1(e) the line passing 
through the origin and y/§ - the square root of the complex conjugate to e (the 
choice of the square root is not important). Careful analysis of the rearrangement 
on Figure [5] shows that F satisfies the following relation: for every edge e S fi 
orthogonal projections of the values of F at its endpoints on the line 1(e) coincide. 
We denote this common projection by F(e) as it would also be given by the same 
formula (fl9|) with z taken on the edge e (to be exact one has to divide by 2 cos(7r/8) 
to arrive at the same normalization). 

It turns out to be a form of discrete analyticity, and implies (but does not 
follow from) the common definition. The latter asks for the discrete version of the 
Cauchy-Riemann equations di a F — id a F to be satisfied. Namely for every lattice 
square the values of F at four corners (denoted u, v, w, z in the counter-clockwise 
direction) should obey 



Remark 13. In the complex plane holomorphic (i.e. having a complex derivative) 
and analytic (i.e. admitting a power series expansion) functions are the same, so 
the terms are often interchanged. Though the term discrete analytic is in wide use, 
in discrete setting there are no power expansions, so it would be more appropriate 
to speak of discrete holomorphic (or discrete regular) functions. 

As discussed above, F solves a discrete version of the Riemann Boundary Value 
Problem (|17|) with a = 1 — Ak = 1/2, which was solved in the continuum case by 



V$'- It remains to show that as the lattice step goes to zero, properly normalized 

F converges to the latter. 

A logical thing to do is to integrate F 2 to retrieve $. Unfortunately, the square 
of a discrete analytic function is no longer discrete analytic and so cannot be 
integrated. However it turns out that there is a unique function H = Im J F 2 dz, 
which is defined on the dual lattice by 



where edge e separates the centers of two adjacent squares, black b and white w. 
After writing (|20|) . one checks that 

1. H is well defined and unique up to an additive constant, 

2. H restricted to white (black) squares is super (sub) harmonic, 

3. H = 1 on (counterclockwise) boundary arc ba and H = on (counterclock- 
wise) boundary arc a6, 

4. The (local) difference between H restricted to white and to black squares 
tends uniformly to zero. 



F(z) - F(v) = i (F(w) - F(u)) . 




H(b)-H(w) = \F(e)\\ 



(20) 
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The properties [TJ [5] are consequences of discrete analyticity: [T] a rather direct 
one, while [2] follows from the identity 



where u is a center of white (black) square with two opposite corner vertices x 
and y (particular choice is unimportant). Definition of F implies the property [3J 
The property [4] easily follows from a priori estimates (namely Kaufman-Onsager- 
Yang results [TU [47]). In principle it should also directly follow from the discrete 
analyticity of F and the property [3l 

We immediately infer that H converges to Im$, and after differentiating and 
taking a square root we obtain the following 

Proposition 14. Suppose that the lattice mesh tj goes to zero and a lattice domain 
flj with boundary points dj, bj converges (in a weak sense, e.g. in Caratheodory 
metric) to a domain Q with boundary points a, b as j —* oo. Then away from the 
boundary there is a uniform convergence: 



Since by Proposition[TT]thc function on the right is a martingale for SLE (16/3), 
convergence of the interface to the Schramm-Lowner Evolution with k = 16/3 
follows. 



6. Conclusion 

At the moment the approach discussed above works only for a (finite) number of 
models. Another notable case when it works is the usual spin representation of 
the Ising model at critical temperature on the square lattice, pictured on Figure [21 
where considering a similar observable (partition function with +1 monodromy, cf. 
(fl9)) ) leads to the martingale (fi"6")) and to Schramm-Loewner Evolution with n = 3. 
Interestingly, exactly the same definition of discrete analyticity arises. 

Analogously, examination of partition function with +1 monodromy at z for 
hexagonal loop models (for all values of n at criticality) suggests its convergence to 
conformal martingale Ijlfip. These considerations lead to a new explanation of the 
Nienhuis' Conjecture[2]for the critical value of x. In this case we firmly believe that 
our method works all the way for n — 1 constructing conformally invariant scaling 
limits for the 0(1) model, but convergence estimates still have to be verified. 

Two parallel methods, with observables related to (fT~5^ an d (fT6]) . seem specially 
adapted to the square lattice and the hexagonal lattice correspondingly. However, 
the main arguments work for a large family of four- and trivalent graphs corre- 
spondingly. So we advance towards establishing the universality conjectures. 

Though only for a few models the conformal invariance was proved, the only 
essential missing step for the remaining ones is discrete analyticity, and it can be 
attacked in a great number of ways. 



AH(u) = ±\F(x)-F(y)\ 2 




bj,z) =1 y/$'{Sl,a,b,z) 
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So from our point of view, the perspectives for establishing conformal invari- 
ance of classical 2D lattice models are quite encouraging. Moreover, we can start 
discussing reasons for universality, and try to construct the full loop ensemble start- 
ing from the discrete picture. The approach discussed above is rigorous, but what 
makes it (and the whole SLE subject) even more interesting is that while borrowing 
some intuition from physics, it gives a new way to approach these phenomena. 
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